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INTRODUCTION 


Under  separate  Prime  Contract  (DA-49-082-OSA-3198) 
IITRI  has  undertaken  a  program  for  development  of  large 
charges  for  use  at  sea  as  seismic  sources.  The  program 
includes  the  development  of  slurry  type  explosives  and 
appropriate  containers.  The  first  at-sea  tests  of  10 
ton  shots  were  performed  off  the  coast  of  Southern 
California  on  February  14  and  21,  1968.  The  explosive 
bodies  were  launched  from  a  ship,  with  delayed  sinking 
to  permit  the  vessel  to  reach  a  safe  stand-off  distance. 
The  purpose  of  these  tests  was  to  obtain  design  infor¬ 
mation  for  use  in  the  development  of  charges  in  the 
fractional  kiloton  range.  Underwater  Systems,  Inc., 
under  Prime  Contrar  L  NOnr  4026(00),  furnished  instru¬ 
mentation  and  manpower  to  monitor  a  number  of  parameters. 
These  include: 

(1)  Yield 

(2)  Detonation  depth 

(3)  Terminal  velocity  of  vehicle 

(4)  Detonation  time 

(5)  Vehicle  attitude 

The  results  obtained  are  discussed  in  the  body  of 
this  report  and  are  summarized  in  Table  I.  Details  of 
the  measurement  procedure  are  given  in  Appendices  A  and  B. 
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Table  I  -  Summary  of  Results 

Yield 

Energy  basis  -  12  to  13  tons  TNT  equivalent 

Bubble  Pulse  Period  basis  -  21  tons  TNT  equivalent 


Shot  1 

Shot  2 

Detonation  Depth  (ft) 

1,825 

3,025 

Terminal  Velocity  { ft/sec 

l  18.9 

18.8 

Detonation  Time  (GMT) 

16  Feb  1968 

21  Feb  1968 

18:01:52.99  +.01 

23:32:08.35  +.01 

Vehicle  Attitude 

Essentially  vertical 

Essentially  vertical 
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RESULTS 


1.  TNT  Equivalence . 

The  yield  of  the  explosive  under  test  is  based  on  a 
comparison  with  TNT.  The  standard  shock  wave  and  bubble  pulse 
equations  given  by  Weston,  Reference  (1) ,  are  used.  A  number 
of  different  measures  of  equivalence  are  available,  depending 
upon  the  portion  of  the  pressure  time  history  which  is  used 
for  the  computation.  Figure  (1)  shows  the  shock  wave  for  shot 
No.  2.  Figure  (2)  shows  the  shock  wave  and  bubble  pulses. 

1.1  Peak  Amplitude  of  the  Shock  Wave .  The  peak  amplitude  of 
the  shock  wave  provides  a  significant  measure  of  TNT  equiva¬ 
lence.  The  peak  level  is  given  by: 


P 

o 


2.16  x  104 


1.13 

psi 


(1) 


W  =  yield  in  pounds 
r  =  range  in  feet 

Since  the  peak  amplitude  of  the  shock  wave  is  independent  of 
detonation  depth,  the  experimentally  measured  levels  of  the 
two  detonations  were  plotted  together  in  Figure  (3).  The 
smooth  curves  for  10  and  13  tons  of  TNT  were  computed  from 
Equation  (1)  . 
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The  peak  signal  varies  as  the  cube  root  of  yield,  and 
a  large  change  in  yield  is  required  to  produce  a  small  change 
in  level.  For  example,  a  10%  change  in  yield  introduces  a 
3 >5%  change  in  signal  level.  Since  the  shock  wave  duration 
is  very  short,  some  rounding  of  the  peak  occurs  due  to  the 
limited  system  bandwidth  of  approximately  zero  to  2,000  Hertz. 
This  tends  to  make  the  measured  levels  somewhat  lower  than 
the  true  levels.  Additionally,  system  noise  introduces  some 
fluctuations  which  can  either  increase  or  decrease  the  measured 
levels.  Both  of  these  effects  can  be  observed  in  Figure  (1). 
The  average  value  of  the  data  points  corresponds  to  a  yield 
of  about  11.5  tons  TNT  equivalent.  If  the  one  low  reading 
for  shot  (2)  at  a  range  of  2,190  feet  is  discounted,  an 
average  yield  of  about  13  tons  TNT  equivalent  is  obtained. 

This  is  our  best  estimate  for  this  measure  of  equivalence. 

A  statistical  analysis  is  unwarranted  since  only  a  few  data 
points  are  available. 

For  a  TNT  equivalence  factor  of  1.3,  the  constant 
2.16  in  Equation  (1)  becomes  2.38,  for  the  actual  weight  of 
the  slurry  explosive  used  by  IITRI . 

1.2  Time  Constant  of  the  Shock  Wave .  A  second  measure  of  the 
TNT  equivalence  is  the  decay  constant  of  the  shock  wave.  The 
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pressure  time  history  of  the  shock  wave  is  represented  by: 

p  =  PQe  t//To  (2) 

where:  Tq  is  the  decay  constant,  and  for  TNT 

is  given  by: 

_  co  ,,1/3  /,,l/3  ,  *-0.22 

Tq  =  58  W  '  (W  /r)  usees  (3) 

The  time  constant  is  also  insensitive  to  yield.  A  10%  change 
in  yield  corresponds  to  a  2*s%  change  in  the  time  constant. 

The  exponential  representation  of  the  signal  is  approximate 
and  holds  to  only  one  or  two  time  constants.  This  measure 
is  therefor  primarily  useful  for  validating  the  TNT  equivalence 
determined  from  the  peak  amplitude  measurement.  By  fitting 
the  theoretical  decay  curve  to  the  observed  signal,  significant 
errors  in  peak  level  due  to  rounding  and/or  noise  can  be  deter¬ 
mined,  and  a  correction  applied.  This  is  illustrated  by 
Figures  (4)  and  (5) .  The  solid  curves  are  tracings  of  the 
average  experimental  signal  with  noise  effects  removed  by  eye 
averaging.  The  dashed  curves  are  based  on  Equations  (2)  and 
(3)  for  a  13  ton  TNT  equivalent  yield,  with  the  level  adjusted 
for  best  fit  to  the  experimental  data  at  t  =  T  ,  and  are  shown 
only  when  they  depart  from  the  experimental  data.  The  experi¬ 
mental  peak  level  is  significantly  lower  than  the  empirical 
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level  for  the  data  of  Figure  (4)  ,  and  substantially  the  same 
for  the  data  of  Figure  (5) .  Significant  corrections  were 
found  for  the  two  data  points  of  Figure  (3)  which  showed  the 
largest  deviation  from  the  theoretical  predictions,  as 
indicated  by  the  arrow  heads.  While  the  absolute  value  of 
these  corrections  are  dependent  upon  the  subjective  process 
of  curve  fitting,  it  is  clear  that  in  both  cases  the  deviation 
from  the  theoretical  prediction  is  reduced.  We  therefore 
conclude  that  this  measurement  of  TNT  equivalence  supports 
the  conclusions  reached  from  the  peak  signal  level. 

1.3  Bubble  Pulse  Period .  The  first  bubble  pulse  period  is 
also  a  useful  measure  of  TNT  equivalence.  When  explosives 
are  compared  to  TNT  it  is  usually  found  that  this  measure 
gives  different  results  than  that  obtained  from  the  shock 
wave  parameters.  The  shock  wave  equivalence  determines  the 
total  available  radiated  energy  at  the  high  frequencies,  and 
a  portion  of  the  energy  at  the  low  frequencies.  The  bubble 
pulse  period  determines  the  spectral  interference  pattern  between 
the  shock  wave  and  bubble  pulse.  A  high  TNT  equivalence  lowers 
the  frequency  at  which  the  first  spectral  maximum  occurs, 
increasing  effectiveness  as  a  seismic  source. 

The  bubble  pulse  period  for  TNT  is  given  by: 

t  =  4.36  W1/3do_5/6  (4) 

where:  dQ  is  the  detonation  depth  plus  32  feet 
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For  the  first  detonation  the  measured  bubble  pulse  period  was 
.287  secs,  and  the  detonation  depth  was  1,825  feet.  For  the 
second  detonation  the  comparable  figures  were  .188  secs  and 
3,025  feet.  From  Equation  (4)  this  yields  a  TNT  equivalence 
of  21  tons  for  both  shots. 

For  a  TNT  yield  equivalence  factor  of  2.1,  the  constant 
4.36  in  Equation  (4)  becomes  5.59,  when  the  actual  weight  of 
the  slurry  explosive  is  used  for  computation. 

1.4  Bubble  Pulse  Parameters .  Additional  measures  of  TNT 
equivalence  are  the  peak  amplitude  and  decay  constant  of  the 
first  bubble  pulse.  The  pressure  time  history  is  represented 
by  an  increasing  exponential  followed  by  a  decreasing  exponential, 
with  identical  time  constants.  These  formulae  are  less  exact 
than  those  used  for  the  shock  wave.  In  particular,  the  pre¬ 
dicted  cusp  is  rarely  observed,  and  the  time  constants  for 
each  side  of  the  signal  are  usually  somewhat  different. 

Equivalence  calculations  are  therefore  less  reliable.  Of 
perhaps  greater  interest  is  the  energy  density  at  the  low 
frequencies.  This  can  be  computed  more  accurately  by  extrapolat¬ 
ing  the  bubble  pulse  peak  to  form  an  exponential  cusp  prior  to 
reading  the  experimental  peak  value  and  time  constant.  When 
this  is  done  we  obtain  a  TNT  equivalent  of  11  tons  based  on 
peak  signal,  and  13  tons  based  on  the  time  constant.  The 
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governing  equations  are : 


=  3,450  (W1/3  /r)  psi 


(5) 


Tx  =  1.39  W 


1.5  Energy.  The  effectiveness  of  the  explosive  as  a  seismic 
source  is  determined  by  the  energy  density  at  low  frequencies. 
For  the  shock  wave  and  first  bubble  pulse,  this  is  given  by: 


2  2 

8pi  V 


pc 


(7) 


(8) 


From  the  TNT  equivalence  determination  made  for  PqI  Tq,  P^, 
and  T^  we  obtain  about  13  tons  TNT  equivalent  for  the  shock 
wave  and  12  tons  for  the  bubble  pulse.  We  therefore  conclude 
that  the  slurry  explosive  used  provides  20%  to  30%  more  energy 
than  TNT.  In  addition,  the  increase  in  the  bubble  pulse  period 
corresponding  to  an  equivalent  110%  increase  in  yield,  decreases 
the  frequency  for  the  first  spectral  maximum.  As  a  result,  the 
effective  energy  available  for  seismic  detection  at  long  range 
can  be  considerably  greater  than  that  for  an  equal  weight  of 
TNT,  depending  on  the  detonation  depth. 
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2.0  Dynamic  Characteristics .  Data  on  sink  rate,  detonation 
time  and  angle  to  the  vertical  were  determined  from  data 
obtained  with  the  acoustic-radio  telemetry  system. 

2.1  Sinking  Rate  and  Detonation  Depth .  Figures  (6)  and 
(7)  show  the  depth  vs.  time  for  the  1st  and  2nd  shots, 
respectively . 

Shot  1. 

Detonation  depth  -  1,825  ft. 

Terminal  velocity  -  18.9  ft/sec 

Shot  2. 

Detonation  depth  -  3,025  ft. 

Terminal  velocity  -  18.8  ft/sec 

During  the  time  periods  when  no  data  is  shown,  the  tele¬ 
metry  signal  was  not  received.  Possible  causes  of  this 
phenomena  are  discussed  in  Appendix  C.  It  is  concluded 
that  the  probable  cause  was  air  bubble  screening  due  to 
venting  around  the  pinger  body  after  the  rupture  diaphragm 
burst  during  descent. 
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APPENDIX  A 


DATA  GATHERING  SYSTEM 


Figure  A-l  shows  the  general  arrangement  of  buoys, 
pingers,  and  recording  vessel  during  the  operation  at 
sea.  Three  data  acquisition  systems  were  used  as  described 
below.  The  acoustic-radio  telemetry  link  was  the  primary 
system. 

1.  Acoustic-Radio  Telemetry 

An  acoustic  pinger  was  attached  to  the  explosive 
vehicle.  Each  of  two  buoys  received  the  pinger  signal  on 
a  hydrophone  and  relayed  the  data  by  radio  telemetry  to 
the  receiving  ship. 

2.  Hard  Wire  -  Radio  Telemetry 

A  buoy  connected  by  torpedo  wire  to  the  pinger 
data  encoding  circuits  provided  a  hard  wire  link,  by¬ 
passing  the  acoustic  telemetry.  The  data  is  relayed  to 
the  recording  vessel  by  radio  telemetry. 

3.  Shipboard  System 

Hydrophones  and  shock  wave  gages  were  used  over 
the  side  at  the  recording  vessel  for  additional  back  up. 
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Acoustic-Radio  Telemetry 


An  acoustic  pinger  operating  at  4.166  kilohertz  was 
mounted  on  each  of  the  explosive  vehicles.  The  source 
level  was  approximately  75  db  re  one  microbar  at  one  yard. 

The  output  consisted  of  a  series  of  coded  pulses. 

a.  An  internal  chronometer  controlled  the  emission 
of  15  millisecond  pulses  at  one  second  intervals. 

b.  In  the  zero  to  1/2  second  time  interval  following 
each  of  the  timing  pulses  an  8  millisecond  pulse  controlled 
by  an  internal  depth  sensor  was  generated.  The  exact  time 
delay  is  proportional  to  vehicle  depth. 

c.  In  the  1/2  to  one  second  time  interval  following 
each  of  the  timing  pulses  a  second  8  millisecond  pulse 
controlled  by  an  internal  yaw  indicator  was  generated. 

The  exact  time  delay  is  related  to  yaw,  as  discussed  in 
Appendix  B. 

Each  of  the  two  acoustic-radio  telemetry  buoys  contained 
a  hydrophone  at  a  depth  of  about  50  feet,  and  a  deep  tour¬ 
maline  gage  to  monitor  the  explosive  signal.  The  two 
buoys  were  launched  prior  to  the  vehicle  launch.  One  buoy 
was  tied  off  from  the  vehicle  with  300  feet  of  line.  A 
release  mechanism  was  inserted  in  the  line  and  is  actuated 
when  the  vehicle  sinks.  The  second  buoy  was  tied  off  from 
the  first  with  200  feet  of  line. 
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The  acoustic  telemetry  signal  is  relayed  to  the 
monitoring  vessel  for  recording.  The  controlled  pulse 
system  permits  a  direct  measurement  of  vehicle  depth  and 
attitude  as  a  function  of  time.  A  typical  data  sequence 
is  shown  in  Figure  A-2. 

Prior  to  sinking,  the  coded  pulse  sequence  controlling 
the  pinger  is  relayed  to  the  recording  vessel  by  radio 
telemetry  from  the  wire-link  buoy  (see  below) .  The 
output  of  a  shipboard  chronometer  and  WWV  time  signals 
are  also  recorded.  Thus,  the  pinger  time  series  can  be 
related  to  absolute  time.  The  shock  wave  resulting  from 
the  detonation  overloads  the  hydrophone  system,  but  the 
arrival  time  is  determined.  From  the  arrival  time  of  the 
preceding  timing  pulse  emitted  by  the  pinger  the  acoustic 
travel  time  is  known,  and  the  detonation  instant  is  det¬ 
ermined.  (See  Appendix  D  for  shock  wave  velocity  correction) . 
From  the  arrival  time  of  the  shock  wave  at  the  tourmaline 
gage  the  range  to  the  gage  is  determined. 
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Hard  Wire  -  Radio  Telemetry 


The  hard  wire  radio  telemetry  buoy  was  mounted  on 
and  launched  with  the  vehicle.  Upon  sinking  the  buoy  is 
submerged  to  a  shallow  depth  until  the  release  mechanism 
actuates  at  which  time  the  buoy  pops  to  the  surface.  One 
reel  of  wire  is  attached  to  the  buoy,  and  a  second  to  the 
vehicle.  As  the  vehicle  sinks  wire  is  payed  out.  The 
wire  connects  the  radio  telemetry  in  the  buoy  to  the  pulse 
control  circuits  in  the  pinger.  Data  on  vehicle  depth 
and  yaw  is  relayed  by  radio  telemetry,  bypassing  the 
acoustic  link.  The  buoy  does  not  have  any  provision  for 
monitoring  the  explosive  signal. 

This  buoy  system  operated  properly  in  harbor  trials. 
During  the  field  operation,  the  release  mechanism  was 
actuated  during  launch  of  the  first  vehicle.  The  buoy 
rolled  off  and  became  inoperative.  This  malfunction  was 
eliminated  by  modification  to  the  release  system  between 
trials.  In  the  second  test  the  buoy  behaved  properly  during 
launch  and  telemetered  data  up  to  the  time  the  vehicle  sank 
after  which  telemetry  signal  was  lost.  Examination  of  the 
buoy  upon  completion  of  the  tests  indicated  that  failure 
was  due  to  water  leakage  at  the  antenna  during  submergence. 
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Shipboard  System 

An  over-the-side  hydrophone  and  two  tourmaline  gages 
for  detecting  the  explosive  signals  were  placed  over  the 
side.  The  primary  purpose  of  this  system  is  to  obtain 
additional  data  on  the  pressure  time  history  of  the  explo¬ 
sive  signal.  The  peak  level  of  the  shock  wave  for  shot  (1) , 
at  a  range  of  approximately  12,000  feet,  shown  in  Figure  (3), 
was  obtained  with  the  shipboard  system. 

Redundancy 

As  is  apparent  from  the  preceding  discussion  that  a 
considerable  amount  of  redundancy  is  built  into  these 
systems.  This  has  been  done  in  order  to  insure  that  sat¬ 
isfactory  data  is  obtained  despite  failure  of  individual 
components  or  entire  sub-systems.  This  redundancy  is 
essential,  since  each  test  is  a  one  shot  operation. 
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IS  60  rf 


Time  (msecs) 


T^glcal  P^lse  time  sequence  as  received 
and  recorded  after  acoustic  and  radio 
telemetry.  This  was  the  last  pulse 

si?eeSehiMCeiVed  fr°m  the  first  explo¬ 
re  vehicle  prior  to  detonation. 
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APPENDIX  B 
YAW  MEASUREMENTS 


To  measure  gross  yaw,  a  tilt  indicator  was  pro¬ 
vided  by  IITRI  and  placed  in  the  pinger.  The  tilt 
indicator  controlled  the  time  delay  between  the  timed 
one  second  acoustic  pulses  and  the  yaw  indicator  pulses. 
The  time  delay  as  a  function  of  tilt  angle  is  shown  in 
Figure  B-l.  As  can  be  noted,  the  linear  portion  of  the 
curve  lies  between  approximately  30°  and  150°.  Since 
IITRI  was  primarily  interested  in  gross  yaw  the  decreased 
sensitivity  between  0°  and  30°  was  acceptable.  Figures 
B-2  and  B-3  show  the  results  obtained  for  the  first  and 
second  vehicles  respectively.  In  both  cases  the  measured 
time  delays  indicate  a  near  vertical  attitude.  The  fluc¬ 
tuations  indicate  reading  accuracy. 
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APPENDIX  C 


LOSS  OF  ACOUSTIC  TELEMETRY  SIGNAL 


During  the  time  period  when  data  is  not  shown,  tele¬ 
metry  signal  was  not  received.  The  cause  is  not  entirely 
understood.  The  several  possibilities  are  discussed  below. 

1.  Ray  path  calculations  based  on  the  sound  velocity 
profile  and  the  pinger-hydrophone  configuration  indicate 
that  good  signal  should  have  been  received.  Acoustic 
shadowing  due  to  refraction  did  not  occur. 

2.  Interference  due  to  direct  and  surface  reflected 
signals  would  introduce  ping-to-ping  fluctuations.  Since 
this  was  not  observed,  and  the  loss  of  signal  occurred  over 
the  same  vehicle  depth  range  for  both  buoys  and  for  both 
shots ,  surface  phenomena  have  been  discounted  as  a  possible 
cause . 

3.  The  acoustic  pinger  was  mounted  on  the  explosive 
vehicle  as  shown  in  Figure  C-l.  Computation  of  the  re¬ 
flected  signals  from  the  vehicle  end  plate  indicates  that 

\. 

interference  between  a  direct  signal  and  a  vehicle  reflected 
signal  would  not  have  caused  loss  of  signal  during  the  in¬ 
terval  in  which  it  occurred. 

4.  Acoustic  shadowing  by  the  vehicle  could  occur  in 
the  event  of  gross  yaw.  Since  the  test  results  indicate 
the  absence  of  gross  yaw  this  possibility  has  been  dis¬ 
counted. 
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5.  It  is  possible  that  loss  of  signal  was  due  to 
air  bubble  screening  of  the  pinger,  but  it  is  not  possible 
to  prove  this  conclusively.  The  pinger  extends  through 
the  vehicle  end  plate  into  the  buoyancy  chamber.  The 
pinger  body  was  bolted  to  the  end  plate.  The  clearance 
space  was  filled  with  RTV  rubber  to  prevent  water  entry 
through  this  space  while  the  vehicle  was  on  the  surface. 

Since  this  was  a  simple  mastic  fill,  without  back-up 
rings,  the  RTV  could  be  squeezed  out  by  a  pressure  head. 

When  the  vehicle  leaves  the  surface  and  sinks,  the  ex¬ 
ternal  water  pressure  will  be  greater  than  the  internal 
pressure  in  the  buoyancy  chamber  until  the  rupture  dia¬ 
phragm  breaks.  During  this  time  interval  water  will  enter 
the  chamber  through  the  vent  holes,  and  if  the  pressure 
differential  becomes  sufficiently  large,  the  RTV  would  be 
squeezed  out.  When  the  rupture  diaphragm  breaks,  the 
rapid  entry  of  water  equalizes  the  pressure  head,  and 
the  remaining  air  vents  under  gravity  through  the  vent 
holes,  and  the  clearance  space  around  the  pinger  body. 

This  would  produce  a  cylindrical  air  bubble  screen  around 
the  pinger  head,  greatly  attenuating  the  radiated  signal, 
until  all  air  had  vented.  This  possible  cause  is  con¬ 
sistent  with  the  observed  loss  of  data  on  both  hydrophones 
for  the  same  range  of  water  depths  for  both  shots,  but  cannot 
be  conclusively  demonstrated. 
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APPENDIX  D 

SHOCK  WAVE  VELOCITY  CORRECTION 

All  range  and  timing  computations  have  been  based 
on  acoustic  propagation  velocities.  Due  to  the  hiah 
pressure  associated  with  the  shock  wave  it  propaqates 
at  a  higher  velocity.  The  measured  detonation  time  can 
be  corrected  for  this  second  order  effect. 

For  the  slurry  explosive  used  in  these  tests  we  have 
found  that  the  peak  pressure  as  a  function  of  the  actual 
weight  of  explosive  used  is  given  by: 

/  1/3  \  1 ' 1 3 

Pq  =  2.38  x  104  — j  ° 


The  propagation  velocity  increases  at  the  rate  of  .04 
ft/sec/psi  for  sea  water.  The  propagation  velocity  of 
the  shock  wave  is  therefore  given  by: 
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The  travel  time  to  a  range  R,  for  propagation  at 


the  acoustic  velocity  is  given  by: 


o 

ior  the  shock  wave,  we  obtain: 
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dr 
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If  we  let  x  equal  the  correction  to  the  shock  wave 
propagation  time  we  obtain: 
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and : 
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Due  to  the  exponent  of  1.13  this  cannot  be  easily 
integrated.  If  we  replace  the  exponent  1.13  by  unity  we 
can  easily  determine  an  upper  limit  for  the  correction 
time  x. 
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c  r  +  954  W 
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which  yields: 
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For  10  tons  of  slurry  at  a  range  of  2,000  feet,  we  obtain: 


x  <  .009  seconds 
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This  correction  has  been  applied  to  obtain  the 
detonation  times  given  in  Table  1. 
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